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We investigate the locality property of the overlap-Dirac operator on gauge configurations gen- 
erated with extra Wilson fermions. By such extra terms we expect that the structure of the Aoki 
phase would change drastically. In particular, we study the possibility of defining the overlap- 
Dirac operator in the strong coupling regime keeping its exponential locality. 
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1. Introduction 

Locality of the overlap-Dirac operator is not obvious, since its definition [p]] 

aD 0V = \ + y 5 -^- (1.1) 

includes an operator Hw in the denominator. Here, Hw is the hermitian Wilson-Dirac operator 
Hw = Y^Dyv, that is used as a kernel to construct D ov . Superficially, if the eigenvalue spectrum of 
H w contains near-zero modes, the overlap-Dirac operator may violate locality. 

The locality is known to be satisfied at weak couplings. To be specific, on background gauge 
configurations satisfying some smoothness condition, one can show that D ov is exponentially local- 
ized [||], i.e. |(D ov )xy| < exp(— |x — y\/£) with a localization length I. This condition is however too 
strong for practical setup used in present lattice QCD simulations; numerical tests are necessary 
for more realistic cases. 

Golterman and Shamir conjectured that the overlap-Dirac operator defined outside of the Aoki 
phase is local [J3|] . Contrary to the original argument [Q], the Aoki phase in this case defined by the 
profile of near-zero modes: inside the Aoki phase the near-zero modes are extended in space, while 
they are localized outside. The value of eigenvalue X c above which the eigenmodes are extended 
is called the mobility edge borrowing the terminology of condensed matter physics. Then, the 
localization length of the overlap-Dirac operator is determined either by 1/A C or l(X) (|A| < A c ), 
where 1{X) is the localization length of the individual low-lying modes. Thus, the question of the 
locality crucially depends on the background gauge field. 

The origin of the near-zero modes of \Hw\ is the roughness of the gauge configuration. A 
simple analytic example is given in [Sj. Therefore, the localization length is expected to increase 
toward strong couplings or coarse lattices, and the definition of the overlap-Dirac operator becomes 
more difficult. To avoid this problem one may introduce additional terms to the lattice action, such 
as those proposed in They consist of two flavors of heavy Wilson fermions and their associated 
ghosts carrying a twisted mass term 

Sex =£zW D w (m ) X (x) + £0 (x) [D w (m ) + ijUy 5 T 3 ] (x) , (1.2) 

X X 

where % denotes Wilson fermions with a negative mass mo. The second term represents bosonic 
fields to cancel the bulk of the effects of Wilson fermions. In fact, the action generates a suppression 
factor 



det 



(1.3) 



H w (m ) 2 + pi 2 _ 

in the partition function. Then, the gauge configuration with small eigenvalues of Hw(mo) lower 
than jit is suppressed, and the potentially dangerous near-zero modes disappear from the eigen- 
value spectrum. Since the low-lying eigenvalue can never cross zero due to the suppression factor 
detH\v(mo) 2 , global topology of the gauge field configuration does not change under continuous 
deformations. The simulations are thus confined in a given topological sector with this action. 

In this work, we investigate the spectrum of Hw and the locality of D ov in the strong coupling 
regime with or without the extra Wilson fermion terms. The purpose of the study is to explore the 
possibility to use the overlap fermion at coarser lattices than currently available ones. 
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A* 





0.2 


P 


5.95 5.83 5.63 5.50 5.43 


5.78 5.68 5.48 5.28 


a(fm) 


0.10 0.12 0.15 0.27 


0.10 0.12 0.14 0.27 



Table 1: Lattices parameters 




40 60 

confid 



P=5.95,|i=0.0 
B=5.83,n=0.0 
0=5.63,11=0.0 
P=5.50,n=0.0 
0=5.43,(1=0.0 
0=5.78,(1=0.2 
0=5.68,(1=0.2 
0=5.48,(1=0.2 
0=5-28,(1=0.2 



Figure 1: Near-zero eigenvalue distribution of \Hw | plotted in the logarithmic scale. The five clusters on the 
left show the results for the standard gauge action while the right four clusters are those with the low-mode 
suppression term. 



2. Eigenvalue distribution 

In this study, we use quenched lattices (no dynamical light fermions) at the j8 values listed 
in Table []]. Lattice size is 16 3 x 32 and the gauge action is the standard Wilson gauge action. We 



consider the lattices with and without the low-mode suppression term (|L3J). In the Table, the lattices 
without that term are denoted by /J. = as the extra factor cancels in this case. The large negative 
mass amo is —1.6, and the number of gauge configurations studied is 10 for each parameter. Lattice 
spacing a is estimated using the Sommer scale ro extracted from the static quark potential. 

We calculate 20-40 lowest eigenvalues of \Hw \ on these lattices. Figure [I] shows those near- 
zero eigenvalue distribution in a logarithmic scale. It is clear that the number of low-modes in- 
creases on coarser (or smaller /3) lattices. With the low-mode suppression term (four right clusters), 
they are indeed highly suppressed. Compared at similar lattice spacings, the lowest eigenvalue is 
1-2 orders of magnitude higher. 

According to the Banks-Casher relation [p||, the absence of the low-lying modes implies that 
the flavor-parity symmetry broken phase as defined by Aoki is not entered for these lattices. This 
is the effect of the low-mode suppressing term, though the gap is rather small at coarser lattices. 
The practical question is then how large the localization length is for these coarse lattices, which is 
addressed in the next section. 
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Figure 2: fhirSf) on the quenched gauge configurations without the low-mode suppressing term (/I = 0.0). 
The data are shown for j3 = 5.95 (top left), /3 = 5.83 (top right), j8 = 5.63 (middle left), j3 = 5.50 (middle 
right), and j3 =5.43 (bottom). The eigenmodes are binned in different ranges of the eigenvalues as indicated 
in the legend of each plot. 

3. Localization of eigenmodes 

In order to investigate the locality of D ov we look at the spatial profile of low-lying modes. For 
each eigenmodes <j>i(x) of Hw, we define p,(x) and fi(r) as 

Pi = $ (*) & (*) , Pi (xo ) = max {pi (x) } , (3.1) 

X 

Mr) = {pi(r)\r=\x-x \} (3.2) 

following to [Q]. Namely p,(x) is the strength of the mode, and fj(r) represents the profile of that 
mode as a function of the distance from the position where p,(x) has a maximum. In calculating 
fi(r), different orientations giving the same r are averaged. 

In Figure Q we plot /i(r) averaged over configurations after binning the eigenmodes in dif- 
ferent ranges of their eigenvalues, that we call fbinif)- At the j8 values above 5.63, we find a clear 
fall-off of the eigenmodes as a function of r. That is true even at the highest bin we measured. This 



4 



Locality of the overlap -Dirac operator 



Yong-Gwi Cho 




- ++ -\ ■v 



























Figure 3: Same as Figure g but with the low-mode suppression term (/x = 0.2). The data are shown for 
J3 = 5.78 (top left), j3 = 5.68 (top right), J3 = 5.48 (bottom left), and j8 = 5.28 (bottom right). 

indicates that D ov constructed on these gauge configurations is local with the length controlled by 
the fall-off of these low-lying modes. 

At /3=5.50 and 5.43, on the other hand, we find that funij) becomes flat beyond r ~ 5. It 
means that the system is already in the Aoki phase. This results is consistent with results of JT6|], 
where the mobility edge falls down to zero at jS = 5.5 and amo = — 1.5. 

Once we introduce the low-mode suppression term, not only the eigenvalue spectrum but the 
profile of the eigenmode changes as shown in Figure |5[ Even at the /3 value as low as 5.28, where 
the lattice spacing is roughly 0.27 fm, the low-lying modes are still localized. This implies that 
these parameter regions are outside of the Aoki phase. In other words, the Aoki phase structure 
is drastically changed by the effect of the low-mode suppressing term. The phase structure of the 
Wilson fermion would thus be changed as illustrated in Figure f|. 

4. Locality of overlap-Dirac operator 

The localization length of the overlap-Dirac operator D ov can also be extracted directly from 
the behavior of D ov . We first set a source field as 

, , f 1 x= (0,0,0,0) 
T] a (x) = { 1 ' \ (4.1) 

[ otherwise 

where a is an index of internal degree of freedom for fermions. Then we calculate a norm of a 
vector y (x) = s g n (Hw) V i x ) 

/(r) = ||^(x)||,r=||x||. (4.2) 
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Figure 4: Expected phase structure of the Wilson fermion after the low-mode suppression term is included. 




Figure 5: (f(r)) as a function of r. Data for the low-mode suppressed gauge configurations: ji = 0.2. 

In Figure [5L we plot (/(r)) as a function of r on the gauge configurations generated with the 
low-mode suppressing term. For all four j8 values calculated, (f(r)) is rapidly decaying as distance 
r increases. 

To determine the localization length, we fit (f(r)) to an exponential function 



</(r))=cexp(-r//). 



(4.3) 



at large distances r. The results for the localization length / is listed in Table ||. It shows that 
the overlap-Dirac operator can be properly defined at strong couplings as far as the low-mode 
suppressing term is introduced. 



p 


/ (localization length) 


5.78 


0.58 


5.68 


0.57 


5.48 


0.60 


5.28 


0.76 



Table 2: Localizationion length calculated on the gauge configurations with the low-mode suppressing term 
/I =0.2. 
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5. Conclusion 

We studied the effect of the low-mode suppressing term on the locality of the overlap-Dirac 
operator. By inspecting the low-lying eigenmodes of Hw, we find the the mobility edge is finite 
even at /? = 5.28, which corresponds to a ~ 0.27 fm. 

Numerical calculations are performed on Hitachi SRI 6000 at High Energy Accelerator Re- 
search Organization (KEK) under a support of its Large Scale Simulation Program (No. 11-05). 
SH is supported in part by the Grant-in-Aid of the Japanese Ministry of Education (No. 21674002). 
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